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1  INTRODUCTION 


In  [1,  4],  the  authors  proposed  and  analyzed  an  interrogation  (inverse  problem)  methodology  based  on  use 
of  an  acoustic  wave  as  a  reflecting  virtual  interface  for  propagating  impulses.  It  is  by  now  well  accepted 
(e.g.,  see  [2,  7,  11,  14])  that  acoustic  pressure  waves  will  interact  with  electromagnetic  signals  in  ways  that 
often  mimic  interfacial  partial  reflcction/partial  transmission  for  the  electromagnetic  waves.  The  response 
of  atomic  electrons  to  an  applied  electrical  field  in  a  material  medium  results  in  a  material  polarization  with 
a  concomitant  index  of  refraction  that  is  a  function  of  the  local  density  in  the  material.  Thus,  material 
density  fluctuations  produced  by  a  sound  wave  induce  perturbations  in  the  index  of  refraction.  Previous 
computational  work  in  [1,  3]  suggested  that  it  might  be  possible  to  detect  reflections  of  microwave  frequency 
EM  waves  from  a  slowly  (relative  to  the  speed  of  the  EM  wave)  moving  acoustic  wave  front.  These  efforts 
focused  on  reflections  in  a  Debye  medium.  The  authors  made  an  argument  for  a  simple  pressure  dependent 
dielectric  model  in  which  the  Debye  parameters  exhibit  a  linear  acoustic  pressure  dependence.  In  [1],  finite- 
element  simulations  for  a  simple  ID  geometry  demonstrated  computationally  that  EM  reflections  from  the 
acoustic  pulse  are  possible.  These  findings  were  confirmed  with  2D  computations  in  [3].  The  results  of 
[1,  3]  consisting  of  a  theoretical  framework  as  well  as  computational  validation  of  such  an  approach  provide 
ample  motivation  for  significant  “proof-of-concept”  experimental  investigations  of  the  proposed  methodology. 
These  prompted  our  preliminary  experiment  on  which  we  report  here  to  look  for  microwave  frequency  EM 
reflections  from  an  acoustic  pulse. 

A  successful  demonstration  of  such  a  technology  is  important.  There  are  many  potential  applications  for 
detection  of  electromagnetic  (EM)  pulses  off  acoustic  waves,  such  as  nondestructive  evaluation  of  materials, 
detection  of  underground  bunkers  and  mines,  and  medical  imaging.  One  generates  the  electrical  property 
distributions  in  the  body  ( Microwave  maps)  with  the  hope  that  such  properties  of  different  bodily  tissues 
are  related  to  their  physiological  state.  With  such  noninvasive  interrogating  techniques  one  can  study 
properties  and  defects  in  biological  tissues  with  very  little  discomfort  to  the  subjects.  Other  potential 
applications  for  such  interrogation  techniques  are  nondestructive  damage  detection  in  structures  where  very 
high  frequency  electromagnetic  pulses  can  be  used  to  detect  the  location  and  width  of  cracks  that  may  be 
present  [5].  Additional  applications  are  found  in  mine,  ordinance  and  camouflage  surveillance,  and  subsurface 
and  atmospheric  environmental  modelling. 

Presently,  ultrasonic  acoustic  waves  are  used  in  medical  imaging  of  the  heart,  fetus,  etc.  Similarly,  EM  waves 
are  employed  in  the  form  of  CT  scans  and  x-rays.  Each  method  has  advantages  that  make  it  the  method  of 
choice  for  visualization  of  particular  body  structures  and  conditions.  One  can  then  imagine  that  combining 
the  two  techniques  and  detecting  EM  reflections  from  acoustic  fronts  might  also  have  its  own  usefulness  for 
particular  scenarios. 

This  report  represents  a  summary  of  our  experimental  investigations  to  detect  the  scattering  of  microwave 
frequency  electromagnetic  waves  by  induced  acoustic  structures.  The  essence  of  this  experiment  consists 
of  launching  electromagnetic  and  acoustic  waves  that  will  be  coincident  in  a  target  and  using  time  domain 
reflectometry  (TDR)  measurements  to  detect  reflections  of  the  electromagnetic  wave  off  the  acoustic  wave 
front. 
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2  EXPERIMENTAL  SET-UP 


The  guiding  wave  structure  is  a  TEM  half  plane  antenna,  as  shown  in  Figure  1.  An  electromagnetic  wave  is 
launched  at  the  feed  (narrow  end)  of  the  antenna  by  means  of  a  fast  rise-time  (200  ps)  input  pulse  from  an 
AVTECH  AV-1030-C  pulse  generator.  The  TEM  antenna  guides  the  spherical  wave  to  a  target  placed  near 
the  aperture  (large  end)  of  the  antenna.  There  are  two  components  to  the  antenna  design:  the  expanding 
and  parallel-plate  portions.  In  the  expanding  portion  of  the  antenna,  the  height  and  width  of  the  antenna 
are  increasing  so  that  the  expanding  spherical  wave  will  more  closely  approximate  a  plane  wave  when  it  hits 
the  target.  The  target  is  placed  in  the  parallel-plate  portion  of  the  antenna,  which  may  be  terminated  with 
resistive  wires  connecting  the  hot  plate  to  the  ground  plane. 


Figure  1:  The  antenna  and  some  of  the  related  equipment.  An  empty  plexiglass  container  is  placed  in  the  parallel  portion  of 
the  antenna  to  illustrate  the  placement  of  the  agar  target.  Resistive  end  termination  wires  can  be  seen  across  the  aperture. 


Figures  2  and  3  are  schematic  diagrams  of  the  antenna.  The  ground  plane  consists  of  a  1/8  inch  thick 
electronic  grade  copper  plate  supported  on  a  wooden  frame,  while  the  hot  plate  is  constructed  from  100 
mesh  copper  wire  cloth  secured  to  a  wooden  frame  with  30  lb.  monofilament  fishing  line.  The  hot  plate  is 

2.75  in.  wide  at  the  feed  of  the  antenna  and  12  in.  wide  at  the  aperture.  The  spacing  between  the  hot  plate 
and  the  ground  plane  is  0.5  in.  at  the  feed  and  12  in.  at  the  aperture. 

The  agar  slab  dimensions  are  13.00  x  11.25  x  1.00  inches,  although  some  measurements  were  performed 
with  a  thicker  (4  in.)  agar  slab.  The  outer  dimensions  of  the  plexiglass  box  containing  the  agar  slab  are 

13.75  x  12.00  x  1.75  inches.  A  one  half  inch  diameter  hole  is  cut  into  the  center  of  the  back  wall  (wide 
end)  of  the  plexiglass  container.  A  hollow  plexiglass  cylinder  (|  inch  length  and  outer  diameter)  is  centered 
on  the  hole  and  projects  out  from  the  container.  The  acoustic  transducer  (ITC-215  kHz  air  transducer)  is 
inserted  into  the  cylinder  and  and  wedged  in  to  make  firm  contact  with  the  back  of  the  agar  slab.  For  the 
measurements  reported  here,  the  agar  slab  was  placed  in  the  antenna  such  that  the  front  of  the  plexiglass 
container  (facing  the  feed)  was  7.25  in.  from  the  aperture  (end)  of  the  antenna.  Therefore,  the  front  surface 
of  the  agar  was  approximately  44.125  in.  from  the  feed  of  the  antenna. 
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TOP  VIEW 


Figure  2:  Top  view  of  antenna  with  agar  slab  (not  to  scale). 
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GROUND  PLANE 


SIDE  VIEW 


Figure  3:  Side  view  of  antenna  with  agar  slab  (not  to  scale). 


A  block  diagram  of  the  experimental  setup  is  shown  in  Figure  4.  Time  domain  reflectometry  (TDR)  measure¬ 
ments  are  used  to  measure  the  reflections  of  the  electromagnetic  pulse  from  structures  inside  the  antenna. 
To  carry  out  the  TDR  measurements  we  used  a  PE2069  50fi  Resistive  Power  Divider  from  Pasternack  Elec¬ 
tronics.  The  power  divider  has  a  frequency  range  of  DC-6  GHz,  an  insertion  loss  of  7  db,  and  a  max  VSWR 
of  1.20:1.  The  50fl  power  divider  splits  the  signal  with  minimal  reflections.  The  input  pulse  from  the  Sig¬ 
nal  Generator  travels  through  an  RG-174  cable  to  the  underside  of  the  ground  plane,  where  it  connects  to 
an  SMA  Female  Bulkhead  connector  (PE4101  Pasternack)  mounted  through  the  ground  plane.  The  inner 
conductor  of  the  connector  feeds  through  the  ground  plane  and  is  connected  to  the  hot  plate  by  means  of  a 
copper  wire  connection.  Referring  to  the  experimental  set-up  in  Figure  4,  it  can  be  seen  that  we  measure 


Figure  4:  Diagram  of  the  TDR  experimental  set-up.  Lo  =  6  ft  and  La  =  6  or  12  ft. 

the  input  pulse  after  the  signal  is  divided  and,  subsequently,  measure  reflected  signals  coming  back  from 
the  antenna.  These  measurements  are  accomplished  using  an  Agilent  Infinium  54854A  oscilloscope  that  is 
triggered  by  a  TTL  signal  from  the  Signal  Generator. 
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2.1  Frequency  Composition  of  Input  Pulse 


A  typical  microwave  input  pulse  to  the  antenna  is  plotted  in  Figure  5.  The  input  pulse  is  approximately 
8  ns  wide  and  has  a  peak  voltage  of  2  V.  Information  about  the  underlying  frequency  composition  of  the 
input  pulse  is  obtained  using  a  fast-fourier  transform  (FFT)  of  the  input  pulse  (Fig.  6).  Plotted  on  a  log-log 
scale,  it  can  be  seen  that  the  input  power  drops  linearly  as  a  function  of  frequency  in  the  range  40  MHz  - 
4  GHz.  Above  4  GHz  there  is  a  steep  drop-off  in  power.  Interestingly,  there  seems  to  be  a  small  burst  of 
power  around  20  GHz. 


Figure  5:  A  typical  input  pulse. 


Figure  6:  FFT  of  the  input  pulse  above. 


3  EMPTY  ANTENNA  RESULTS 

We  first  investigated  reflections  from  structures  in  the  antenna  in  the  absence  of  an  acoustic  target  using 
TDR  measurements.  Figure  7  compares  data  for  the  two  cases  of  La  =  6  and  12  ft  (see  Fig.  4).  The  initial 
pulses  (t  =  0)  for  the  La  =  6  and  12  ft  cases  coincide,  as  expected,  representing  the  portion  of  the  signal 
generator  pulse  that  travels  directly  from  the  splitter  to  the  oscilloscope.  The  (time)  positions  of  subsequent 
pulses  depends  upon  the  coaxial  cable  length  La- 

The  first  reflected  pulse  occurs  at  times  t\  =  19  and  38  ns  for  the  6  ft  and  12  ft  cables,  respectively.  This 
pulse  occurs  because  of  reflections  at  the  coaxial  cable/antenna  connection  and  indicate  imperfect  impedance 
matching  between  the  feed  and  the  50H  coaxial  cable.  The  pulse  reflected  at  the  coax/antenna  junction  has 
the  same  orientation  as  the  incident  pulse.  The  reflection  coefficient  is  0.186  (incident  voltage/reflected 
voltage) . 
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The  second  reflected  pulse  corresponds  to  reflections  of  the  transmitted  pulse  from  the  end  of  the  antenna. 
Reflections  from  the  open-ended  antenna  maintain  the  same  orientation.  After  the  pulse  reflects  from  the  end 
of  the  antenna  a  portion  is  transmitted  through  the  feed  to  the  oscilloscope  (incident  time  t2  =  +  while 

the  rest  is  re-radiated  in  the  antenna,  but  with  an  inverted  portion.  Subsequent  reflected  pulses  alternate  in 
orientation,  each  occurs  at  a  time  tA  after  the  next. 

The  time  of  arrival  of  the  reflected  pulses  can  be  calculated  using  information  about  the  speed  of  propagation 
in  the  RG174  coaxial  cables  (1.9  x  108  m/s)  and  in  air  (3.0  x  108  m/s)  and  the  dimensions  of  the  antenna. 
Table  2  tabulates  calculated  values  for  time  of  arrival  of  reflection  pulses,  where  it  can  be  seen  that  the 
calculated  values  are  in  good  agreement  with  measurements  (Fig.  7). 


052604  data.  PE2069  50  ohm  Tee  data.  1  ave.  No  agar.  No  coax  sensors. 


Time  (ns) 


Figure  7:  TDR  experiment  for  the  two  cases  of  a  6  ft  and  12  ft  cable  from  the  power  divider  to  the  antenna.  Calculated 
reflection  times  are  indicated  by  the  arrows. 


Table  2 .Round-trip  Times 


PULSE  ARRIVAL  TIMES 
(CALCULATED) 

6  ft  COAX  12  ft  COAX 
(ns)  (ns) 


tcOAX 

19.25 

38.5 

tcOAX  +  tA 

27.71 

46.96 

tcOAX  +  2tA 

36.17 

55.42 

tcOAX  +  3  tA 

44.63 

63.88 

tcOAX  +  4  tA 

53.09 

72.34 
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3.1  Radiation  Absorbing  Tiles 


Figure  8  compares  the  reflected  pulses  with  and  without  radiation  absorbers  (Braden  Shielding  Systems)  at 
the  aperture  of  the  antenna.  No  significant  difference  can  be  seen  between  the  two  cases.  The  combination 
of  ferrite  tiles  and  pyramidal  ferrite  absorbers  that  was  used  provides  a  broad  range  of  power  absorption 
in  the  frequency  range  of  the  antenna.  Therefore,  these  results  suggest  that  the  radiated  portion  of  the 
electromagnetic  pulse,  which  is  being  absorbed  by  the  tiles,  is  largely  independent  of  the  reflected  signal  in 
the  TDR  measurements.  Despite  having  little  affect  on  the  reflected  signal,  the  radiation  tiles  were  placed 
at  the  aperture  of  the  antenna  for  all  the  EM/acoustic  interaction  data  reported  here. 


052604  data.  PE2069  50  ohm  Tee  data.  1  ave.  No  agar.  No  coax  sensors. 


Figure  8: 


TDR  measurements  with  (blue)  and  without  (red)  radiation  absorbing  tiles  at  the  aperture. 


3.2  Resistive  End  Termination 


In  the  next  set  of  experiments  we  seek  to  reduce  the  reflections  inside  the  antenna  by  varying  the  end 
termination  at  the  aperture  end  of  the  antenna.  Four  or  five  resistors  (5%  tolerance,  1/4  Watt)  are  connected 
in  parallel  from  the  hot  plate  to  the  ground  plane.  Figure  9  plots  the  TDR  results  for  different  values  of 
(total)  resistance  between  the  hot  plate  and  the  ground  plane.  The  first  peak  in  Fig.  9(a)  is  the  (split)  pulse 
directly  from  the  Signal  Generator.  The  second  peak  is  the  reflection  of  the  other  half  of  the  split  pulse 
at  the  coaxial  cable/antenna  feed  interface.  As  expected,  this  peak  is  independent  of  the  end  termination. 
The  third  peak/ valley  and  subsequent  features  represent  repeated  reflections  of  the  pulse  from  the  end  of 
the  antenna.  These  features  will  be  greatly  affected  by  the  end  termination,  since  reflections  from  the  end 
of  the  antenna  will  be  minimized  when  the  end  termination  resistance  matches  the  characteristic  impedance 
of  the  antenna  at  the  aperture.  From  Fig.  9  it  is  evident  that  the  end  reflections  are  minimized  with  an  end 
termination  of  approximately  125  Cl. 

These  results  suggest  that  the  characteristic  impedance  of  the  antenna  at  the  aperture  end  is  approximately 
125  Cl.  The  antenna  was  designed  to  have  a  377  Cl  characteristic  impedance  at  the  aperture  in  order  to  be 
impedance  matched  to  the  characteristic  impedance  of  radiation  in  free  space.  This  design  was  based  on 
the  assumption  that  the  characteristic  impedance  of  the  end  configuration  most  closely  matched  that  of  a 
parallel  plate  configuration  (Fig.  10a).  Because  the  ground  plane  extends  out  beyond  the  hot  plate,  our 
configuration  may  more  closely  resemble  that  of  a  microstrip  (Fig.  10b).  Using  the  microstrip  formula,  the 
characteristic  impedance  of  the  antenna  at  the  aperture  end  ( H  =  12  and  W  =  12  inches)  is  expected  to  be 
Zq  =  125U,  which  agrees  well  with  the  end  termination  resistance  results  (above). 
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060204  data.  TDR.  1058  ave.  End  termination  comparison. 


Figure  9:  TDR  results  with  different  resistive  end-terminations  at  the  antenna  aperture. 


TOP  FRINGING 

CONDUCTOR  FIELD 


W  ►  GROUND  PLANE 


Z0=r|  H/W 


Z0=r|/(W/H  +  2) 


(a)  Parallel  Plate. 


(b)  Microstrip. 


Figure  10:  Two  different  configurations  that  approximate  the  characteristic  impedance  of  the  antenna  at  the  aperture  end 
The  dimensions  of  the  antenna  are  H  =  12  and  W  =  12  inches  at  the  aperture.  In  the  equations  7/  =  \J //q  / ro  =  377 O  is  the 
characteristic  impedance  of  free  space. 
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ACOUSTIC  WAVES 


4.1  Transducer  Characterization 


The  ultrasonic  transducer  also  acts  as  a  receiver  of  acoustic  waves.  The  receiver  capability  can  be  used  to 
check  that  the  transducer  is  functioning  properly  and  to  measure  the  velocity  of  sound  in  the  agar  target. 
Figure  11  plots  the  signal  that  is  received  by  the  transducer  as  a  function  of  time  for  different  target  distances 
( D ).  The  initial  portion  of  the  signal  is  the  same  for  all  cases  and  represents  the  ring-down  of  the  pulse  in 
the  transducer.  The  secondary  pulses  represent  the  reflection  of  the  acoustic  signal  from  a  hard  target.  As 
expected,  the  locations  (round-trip  times)  of  the  reflected  pulses  increase  as  the  target  distance  increases. 


TRANSDUCER  RECEIVER  CIRCUIT  OUTPUT 


TIME  (ms) 

Figure  11:  Output  of  the  received  transducer  signal  (Volts)  as  a  function  of  time  for  different  target  distances  ( D ). 

The  time  of  arrival  of  the  reflected  pulse  is  plotted  as  a  function  of  the  target  distance  in  Figure  12  (circles) 
along  with  the  calculated  round-trip  time  (line).  Offsetting  the  calculated  times  by  0.11  ms  gives  good 
agreement  between  the  observed  and  calculated  arrival  times.  The  offset  time  probably  represents  a  fixed 
delay  time  in  the  circuit. 


ROUND-TRIP  TIME  VS  TARGET  DISTANCE 


Figure  12:  Observed  (circles)  and  calculated  (line)  arrival  times  for  the  reflected  acoustic  pulses  as  a  function  of  the  target 
distance  D. 


4.2  The  Velocity  of  Acoustic  Waves  in  Agar 


To  measure  the  velocity  of  the  acoustic  waves  in  agar  we  mount  the  transducer  on  the  top  (narrow)  surface 
of  agar  and  point  it  down  into  agar.  The  acoustic  pulse  travels  0.29  m  to  the  bottom  plexiglass  surface, 
where  it  is  reflected  back  to  the  transducer.  The  reflected  signal  received  by  the  transducer  is  shown  in 
Figure  13  (black  curve).  The  initial  portion  of  the  signal  (0-0.5  ms)  represents  the  ring-down  period  of  the 
transducer  receiver  circuit.  There  are  many  oscillations  in  the  signal  from  the  receiver  circuit,  but  distinct 
reflection  “pulses”  are  evident.  It  is  easier  to  look  at  a  min/max  smoothed  plot  of  the  data  shown  above. 


Reflection  of  acoustic  signal  travelling  through  agar 


TIME  (ms) 


Figure  13:  Oscilloscope  measurement  of  transducer  receiver  circuit  signal  for  the  case  of  an  acoustic  wave  travelling  through 
agar.  Reflections  from  the  opposing  plexiglass  surface  can  be  seen  at  0.5  ms  and  at  approximately  0.4  ms  intervals  thereafter. 
The  red  and  green  curves  show  the  min/max  smoothed  curved  for  this  data  ( N  =  100). 

At  each  time  point  the  plotted  value  equals  the  maximum  (red  curve)  or  minimum  (green  curve)  value  of 
the  current  point  and  N  points  ahead  and  behind  the  current  point  (Figure  14  with  N  =  100).  The  velocity 
of  the  acoustic  wave  in  agar  is  calculated  to  be  the  round-trip  distance  in  agar  divided  by  the  time  between 
reflections.  Here  the  velocity  is  measured  to  be  1450  m/s,  which  is  consistent  with  published  values. 
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Reflection  of  acoustic  signal  travelling  through  agar 


Figure  14:  The  min/max  smoothed  plot  of  the  signal  from  the  transducer  receiver  circuit  showing  reflections  of  the  acoustic 
pulse  from  the  plexiglass  surface. 


4.3  Comparison  of  the  Transducer  Receiver  Signal  in  Agar  and  Air 


The  transducer  is  mounted  against  the  back  of  the  agar  as  is  the  case  when  measuring  EM/ Acoustic  in¬ 
teractions.  The  acoustic  pulse  travels  1  inch  before  reflecting  from  the  opposing  plexiglass  surface.  The 
transducer  is  aimed  at  a  distant  point  outside  the  antenna  (in  air).  No  reflections  are  evident  in  the  air, 
because  of  the  large  distance.  The  reflections  in  agar  are  spaced  at  50  /xs  intervals.  The  spacing  should  be 
35  /x s. 


MIN/MAX  SMOOTHED  RECEIVER  SIGNAL  IN  AGAR  SAMPLE  AND  AIR 


Figure  15:  The  signal  from  the  transducer  receiver  circuit  for  a  pulse  traveling  in  air  (blue  and  cyan  curves)  as  compared  to 
a  pulse  traveling  in  the  agar  (red  and  green  curves). 
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5  SYNCHRONIZATION  OF  EM  AND  ACOUSTIC  WAVES 


In  order  to  have  the  electromagnetic  and  acoustic  waves  coincident  in  the  antenna  we  must  have  allow  a 
sufficient  delay  for  launching  the  electromagnetic  wave  after  the  acoustic  wave  is  launched,  since  the  velocity 
of  the  electromagnetic  wave  is  much  greater  than  that  of  the  acoustic  wave. 

The  synchronization  is  accomplished  with  two  timer  circuits  and  a  variable  delay  time  control  on  the  Signal 
Generator  (Fig.  4).  The  first  timing  circuit  triggers  the  second  timing  circuit  as  well  as  providing  the  trigger 
input  for  the  transducer  circuit.  The  second  triggering  circuit  serves  as  an  external  trigger  to  the  Signal 
Generator.  Synchronization  of  the  acoustic  and  EM  pulses  is  achieved  by  selecting  an  appropriate  value  for 
the  delay  time  on  the  Pulse  Generator.  The  delay  time  is  the  time  between  the  arrival  of  the  External  Trigger 
pulse  and  the  output  of  the  EM  pulse  from  the  Pulse  Generator  (Fig.  16).  1  Since  the  speed  of  the  acoustic 
pulse  in  the  agar  is  approximately  17  ^s/inch  and  the  transducer  trigger  pulse  is  134  /is,  the  variable  delay 
should  be  approximately  142  /is  for  the  EM  pulse  to  encounter  the  acoustic  wave  at  the  middle  of  the  1  inch 
thick  agar. 


134  (is 


TRANSDUCER  TRIGGER 


50  |IS 


i-6  (is  CSS  If  < 


TRANSDUCER  PULSE 


PG  EXTERNAL  TRIGGER 


VARIABLE 

< -  DELAY 


EM  PULSE 


Figure  16: 


Schematic  diagram  of  the  signal  timing.  PG  refers  to  the  Pulse  Generator. 


1More  accurately,  the  Delay  Time  is  the  time  between  output  of  the  SYNCH  pulse  and  output  of  the  EM  pulse.  There  is  a 
fixed  propagation  delay  of  60  ns  between  the  External  Trigger  and  SYNCH  pulses. 
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6  EM/ACOUSTIC  REFLECTIONS  IN  AGAR 


In  these  experiments  we  make  TDR  measurements  with  and  without  adding  acoustic  pulses.  Each  waveform 
plotted  is  the  average  of  4096  measurements.  The  waveforms  with  and  without  the  acoustic  pulse  are 
subtracted  to  look  for  differences  that  may  be  due  to  a  reflection  of  the  EM  pulse  from  the  acoustic  wave 
front.  Radiation  absorbing  tiles  and  125  f 1  resistive  termination  were  used  in  the  experiments  reported  in 
this  section. 


6.1  Experiments  With  Different  Delays 


We  next  perform  TDR  experiments  with  an  agar  target  and  acoustic  waves.  We  vary  the  delay  time  to 
sample  a  range  of  possible  positions  of  the  acoustic  wave  in  the  agar  at  the  time  of  incidence  of  the  EM  wave. 
The  speed  of  the  EM  pulse  in  the  agar  is  estimated  to  be  3.4  x  107  m/s,  while  the  speed  of  sound  in  the  agar 
was  measured  to  be  1450  ms/s.  Figure  17  plots  the  positive  portions  of  the  transducer  receiver  signals  for 
the  different  delay  times  that  were  used.  Time  t  =  0  indicates  the  time  at  which  the  EM  wave  is  launched 
in  the  antenna.  As  seen  in  Fig.  17,  the  transducer  receiver  signals  are  shifted  relative  to  the  launch  times 
(t  =  0),  according  to  the  different  delay  times  used. 


MAX-SMOOTHED  (200)  RECEIVER  SIGNAL  IN  AGAR 
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Figure  17:  The  positive  portion  of  the  min/max  smoothed  signal  from  the  transducer  receiver  circuit  for  different  delay  times: 
540  fis  (cyan),  435  (green),  395  (black),  255  (magenta),  174  (red),  and  142  fis  (blue).  The  time  axis  is  measured  relative  to  the 
launch  time  of  the  EM  pulse  (t=0). 
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The  TDR  measurements  for  the  six  different  delay  times,  with  and  without  acoustic  waves,  are  plotted  in 
the  upper  half  of  Figure  18.  Because  the  differences  among  the  measurements  are  small,  all  twelve  curves 
appear  to  be  coincident  at  this  scale.  The  initial  reflection  of  the  EM  pulse  off  the  agar  surface  is  indicated 
by  the  first  negative  swing  in  the  TDR  signals  at  approximately  26.5  ns  (top  plot).  The  first  reflection  of 
the  EM  pulse  from  the  acoustic  wave  is  expected  to  occur  shortly  after  this  time.  Because  of  reflections  from 
other  structures  in  the  antenna,  reflections  of  the  EM  wave  off  the  acoustic  wave  structure,  if  they  exist,  are 
expected  to  be  small  in  comparison.  For  this  reason  we  plot  a  subtracted  signal  in  the  lower  half  of  Figure 
18.  The  subtracted  signal  is  the  difference  between  TDR  measurements  with  and  without  an  acoustic  wave 
in  the  agar  target.  Ideally,  the  plot  of  the  subtracted  signals  would  indicate  only  those  events  related  to  the 
reflection  of  the  EM  pulse  from  the  acoustic  wave.  However,  there  appear  to  be  three  main  features  visible 
in  the  subtracted  signal:  the  large  spikes  associated  with  the  rise  and  fall  of  the  input  pulse,  a  smaller  pair 
of  spikes  that  occur  approximately  19  ns  later,  and  the  presence  of  long-term  oscillations  (60-180  ns). 


091704  data.  TDR.  4096  ave. 
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Figure  18:  The  top  figure  plots  the  TDR  measurements  for  the  different  delay  times  with  and  without  acoustic  signals.  There 
is  little  difference  among  the  measurements,  thus  they  appear  as  a  single  curve.  The  lower  plot  shows  the  subtracted  signals 
(measurement  without  acoustic  pulse  minus  measurement  with  acoustic  pulse)  for  the  different  delay  times:  540  fis  (cyan),  435 
(green),  395  (black),  255  (magenta),  174  (red),  and  142  /xs  (blue).  Note  that  the  units  of  the  lower  plot  is  in  millivolts. 
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Input  pulse  jitter.  In  Figure  19  we  take  a  closer  look  at  the  large  spikes  in  the  subtracted  signal  that  are 
associated  with  the  rise  and  fall  of  the  input  pulse.  Without  jitter,  the  subtracted  voltage  should  be  zero  for 
the  duration  of  the  pulse  width.  The  spikes  in  the  subtracted  signals  that  occur  at  the  leading  and  trailing 
ends  of  the  input  pulse  may  be  due  to  jitter  in  the  pulse  width  or  jitter  in  the  synchronization  signal  that 
triggers  the  oscilloscope.  There  seems  to  be  a  correlation  in  the  negative  and  positive  spikes;  a  negative 
spike  in  the  leading  edge  of  the  pulse  is  generally  followed  by  a  positive  spike  in  the  trailing  edge  of  the  input 
pulse  subtracted  signal.  This  suggest  that  jitter  in  the  synchronization  signal  may  be  the  main  source  of  the 
discrepancies  in  the  subtracted  signal.  If  the  timing  of  the  synchronization  signal  varies  a  little  bit,  the  time 
at  which  the  input  pulses  are  sampled  by  the  oscilloscope  will  vary  a  little  bit.  This  would  produce  spikes 
at  the  edges  of  the  input  pulse  when  the  two  signal  are  subtracted,  since  they  would  be  slightly  offset. 


091604  data.  TDR.  4096  ave. 


Time  (ns) 


Figure  19:  This  figure  gives  an  expanded  view  around  the  time  that  the  input  pulse  arrives.  The  top  figure  plots  the  TDR 
measurements  for  the  different  delay  times  with  and  without  acoustic  signals.  There  is  little  difference  among  the  measurements, 
thus  they  appear  as  a  single  curve.  The  lower  plot  shows  the  subtracted  signals  (measurement  without  acoustic  pulse  minus 
measurement  with  acoustic  pulse)  for  the  different  delay  times:  540  fis  (cyan),  435  (green),  395  (black),  255  (magenta),  174 
(red),  and  142  fis  (blue). 
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Events  in  the  vicinity  of  the  EM/Acoustic  interaction.  Figure  20  also  provides  a  closer  look  at  the 
above  plot  in  the  time  region  where  we  expect  the  EM  and  acoustic  waves  to  be  coincident.  The  speed  of 
the  EM  pulse  in  the  agar  is  3.4  x  107  m/s,  so  the  EM  pulse  should  encounter  the  acoustic  wave  at  the  27-28 
ns  mark.  We  do,  in  fact,  see  voltage  spikes  in  the  subtracted  signal  at  25  ns.  If  the  events  at  the  25  ns  mark 
were  due  to  EM  reflections  from  the  acoustic  wave,  we  would,  because  of  the  different  delay  times,  expect 
these  reflection  events  occur  at  slightly  different  times.  We  can  also  see  from  Fig.  20  that  there  are  voltage 
spikes  in  the  subtracted  signal  that  occur  approximately  6  ns  earlier,  at  the  19  ns  mark.  The  time  at  which 
these  pairs  of  spikes  occur,  19  ns  after  the  leading  and  trailing  edges  of  the  input  pulse,  and  the  temporal 
separation  of  these  spikes  (6  ns)  suggest  that  they  are  probably  related  to  the  reflection  of  the  input  pulse 
at  the  feed  of  the  antenna.  The  variability  in  the  input  pulse  synchronization  is  also  seen  in  its  reflection  at 
the  feed.  We  are  observing  the  input  pulse  jitter  again,  although  it  has  been  attenuated. 


091604  data.  TDR.  4096  ave. 


Figure  20:  This  figure  gives  an  expanded  view  around  the  time  that  the  input  pulse  first  encounters  the  acoustic  wave.  The 
top  figure  plots  the  TDR  measurements  for  the  different  delay  times  with  and  without  acoustic  signals.  There  is  little  difference 
among  the  measurements,  thus  they  appear  as  a  single  curve.  The  lower  plot  shows  the  subtracted  signals  (measurement 
without  acoustic  pulse  minus  measurement  with  acoustic  pulse)  for  the  different  delay  times:  540  fis  (cyan),  435  (green),  395 
(black),  255  (magenta),  174  (red),  and  142  p,s  (blue). 


6.2  Experiment  With  Different  Pulse  Widths 


To  further  explore  the  phenomena  that  we  observed  in  the  above  data,  another  set  of  measurements  was 
performed  where  the  delay  time  was  kept  constant,  but  the  width  of  the  input  pulse  was  varied  (Figure  21). 
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091 704  data.  TDR.  4096  ave. 


Figure  21:  The  top  figure  plots  the  TDR  measurements  for  the  different  input  pulse  widths,  with  and  without  acoustic  signals. 
The  lower  plot  shows  the  subtracted  signals  (measurement  without  acoustic  pulse  minus  measurement  with  acoustic  pulse)  for 
the  different  pulse  widths:  7  (blue),  9  (red),  13  (green),  and  16  ns  (black).  Note  that  the  units  of  the  lower  plot  is  in  millivolts. 

Figure  22  expands  the  view  around  the  input  pulse.  Again  we  can  see  voltage  spikes  in  the  subtracted  view 
(lower  plot)  associated  with  the  leading  and  trailing  edges  of  the  input  pulses.  However,  in  the  case  of  the 
13  ns  pulse  (green  curve),  timing  of  the  acoustic  and  non-acoustic  input  pulses  is  closely  matched  and  the 
subtracted  signal  has  only  very  small  spikes  in  the  vicinity  of  the  leading  and  trailing  edges. 


091704  data.  TDR.  4096  ave. 


Time  (ns) 


Figure  22:  An  expanded  view  of  the  input  pulses  and  the  subtracted  waveforms.  The  top  figure  plots  the  TDR  measurements 
for  the  different  input  pulse  widths,  with  and  without  acoustic  signals.  The  lower  plot  shows  the  subtracted  signals  (measurement 
without  acoustic  pulse  minus  measurement  with  acoustic  pulse)  for  the  different  pulse  widths:  7  (blue),  9  (red),  13  (green),  and 
16  ns  (black).  Note  that  the  units  of  the  lower  plot  is  in  millivolts. 


16 


Figure  23  expands  the  view  in  Fig.  21  in  the  vicinity  of  the  EM/acoustic  wave  reflections.  In  Fig.  18  above, 
we  saw  that  the  events  at  the  25  ns  mark  were  coincident  even  though  we  had  different  delay  times.  In  Fig. 
23  (lower  plot)  we  can  see  that  these  events  are  no  longer  coincident  when  the  pulse  width  is  varied  (constant 
delay  time),  further  suggesting  that  these  events  are  related  to  the  input  pulse  jitter. 


091 704  data.  TDR.  4096  ave. 


Figure  23:  An  expanded  view  in  the  vicinity  of  the  EM/Acoustic  interaction.  The  top  figure  plots  the  TDR  measurements  for 
the  different  input  pulse  widths,  with  and  without  acoustic  signals.  The  lower  plot  shows  the  subtracted  signals  (measurement 
without  acoustic  pulse  minus  measurement  with  acoustic  pulse)  for  the  different  pulse  widths:  7  (blue),  9  (red),  13  (green),  and 
16  ns  (black).  Note  that  the  units  of  the  lower  plot  is  in  millivolts. 


If  we  plot  the  times  at  which  these  events  occur  as  a  function  of  the  pulse  width,  we  see  a  strong  linear 
relationship  (Fig.  24).  In  this  plot,  there  is  no  data  point  corresponding  to  the  13  ns  pulse  width,  since 
there  was  no  clearly  observable  voltage  event  in  the  vicinity  of  the  expected  time  of  the  EM/acoustic  wave 
reflection.  This  is  further  evidence  that  the  observed  voltage  events  are  related  to  the  input  pulse  jitter, 
since  there  was  very  little  jitter  in  the  13  ns  input  pulse  (Fig.  22). 


Figure  24: 


Location  of  the  spikes  in  the  subtracted  signal  (Fig.  23)  as  a  function  of  the  input  pulse  width. 
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6.3  Experiment  With  Thicker  Agar 


The  time  for  the  EM  pulse  to  traverse  the  1  in.  thick  agar  is  roughly  0.75  ns.  This  makes  it  difficult  to 
separate  the  EM  reflection  at  the  front  of  the  agar  from  the  EM  reflection  at  the  acoustic  wavefront.  Figure 
25  plots  the  results  of  measurements  made  with  a  4  inch  thick  agar  target.  As  in  previous  experiments  we 
take  data  with  and  without  an  acoustic  wave  in  the  agar,  as  well  as  recording  a  third  data  set,  also  without 
an  acoustic  signal.  When  plotted  together  all  three  waveforms  appear  indistinguishable  (Fig.  25,  top  plot). 
If  one  of  the  non-acoustic  data  sets  is  subtracted  from  the  acoustic  data  set  we  get  results  that  are  similar  to 
what  we  have  already  seen  (Fig.  25,  bottom  plot,  blue  curve).  Subtracting  one  non-acoustic  data  set  from 
the  other  provides  a  control  for  comparison  to  the  acoustic  data  (Fig.  25,  bottom  plot,  red  curve). 

Both  the  acoustic  and  the  control  wave  forms  in  Fig.  25  (lower  plot)  have  spikes  at  19  and  25  ns,  further 
evidence  that  this  is  a  remnant  of  the  input  pulse  and  not  an  acoustic  reflection  event.  There  seems  to  be 
jitter  in  the  input  pulse  peak  voltage.  Both  sets  of  subtracted  data  are  nonzero  during  the  input  pulse  time. 


1 021 04  data.  TDR.  4096  ave. 


Time  (ns) 


Figure  25:  TDR  and  subtracted  data  with  a  4  inch  thick  agar  target.  In  the  lower  plot,  the  blue  curve  is  the  subtracted  data 
for  acoustic  minus  no-acoustic  measurements.  The  red  curve  is  the  subtracted  data  for  two  no-acoustic  measurements. 
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6.4  EM/Acoustic  Reflections  Without  Agar 


One  feature  of  the  subtracted  TDR  signals  that  is  consistently  observed  is  the  presence  of  long  term  oscilla¬ 
tions.  These  oscillations  are  present  in  the  control  (Fig.  25,  bottom  plot,  red  curve),  although  the  amplitude 
of  the  oscillations  is  much  reduced.  In  order  to  explore  the  possiblity  that  these  oscillations  are  an  artifact 
of  the  experiment  and  not  related  to  EM/acoustic  reflections,  we  measurements,  as  above,  but  without  an 
agar  target  in  the  antenna.  In  this  experiment,  the  transducer  is  mounted  through  a  styrofoam  block  (2.5” 
wide,  10”  high,  and  1”  thick)  that  acts  as  a  support  for  the  transducer.  The  back  of  the  styrofoam  block  is 
flush  with  the  pyramidal  radiation  absorbing  tiles  at  the  aperture  of  the  antenna.  The  front  surface  of  the 
transducer  is  flush  with  the  front  surface  of  the  styrofoam  block,  about  1.5”  from  the  back  of  the  antenna, 
and  5.25”  above  the  ground  plane.  The  transducer  is  pointed  toward  the  feed  of  the  antenna.  The  transducer 
wire  feeds  through  a  small  opening  in  the  panel  of  radiation  absorbing  tiles. 

Figure  26  plots  subtracted  TDR  data  for  the  following  experimental  conditions.  The  subtracted  data  plotted 
with  the  red  curve  was  obtained  in  a  manner  similar  to  the  acoustic  experiments  with  an  agar  target,  that 
is,  TDR  data  with  no  acoustic  wave  is  subtracted  from  TDR  data  taken  with  an  acoustic  wave.  As  is  usually 
done,  the  nonacoustic  data  is  taken  with  the  timing  circuit  disconnected  from  the  transducer  circuit  so  that 
no  trigger  events  are  sent  to  the  transducer  drive  circuit.  The  timing  circuit  is  externally  triggering  the 
signal  generator  (this  data  was  taken  with  a  180  yts  delay  time).  As  the  red  curve  illustrates,  there  are  strong 
long-term  oscillations  similar  to  what  has  been  seen  with  an  agar  target. 

The  data  plotted  with  the  blue  curve  is  the  subtracted  data  of  two  nonacoustic  TDR  measurements,  where 
the  nonacoustic  experimental  conditions  are  achieved  in  two  different  ways.  In  one  case  the  timing  circuit 
is  disconnected  from  transducer  circuit,  as  above.  In  the  second  case  the  transducer  drive  circuit  remains 
connected  to  the  timing  circuit,  but  there  is  no  power  driving  the  timing  circuit  and  the  Signal  Generator 
is  set  at  a  100  Hz  repetition  rate.  As  seen  in  Figure  26,  the  nonacoustic  subtracted  data  shows  similar 
long-term  oscillations  as  the  acoustic/nonacoustic  subtracted  data. 
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Figure  26:  Subtracted  signals  for  TDR  data  taken  with  no  agar  target.  The  blue  curve  is  the  subtracted  signal  from  two 
nonacoustic  data  sets,  taken  under  different  experimental  conditions.  The  red  curve  is  the  subtracted  signal  from  acoustic  and 
nonacoustic  data  sets  without  an  agar  target. 
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For  comparison  purposes,  Figure  27  plots  the  subtracted  data  obtained  under  the  two  nonacoustic  exper¬ 
imental  conditions  along  with  acoustic/nonacoustic  subtracted  data  obtained  with  an  agar  target  in  the 
antenna.  The  data  in  the  lower  plot  of  Figure  27  show  that  the  magnitude  of  the  long  term  oscillations  are 
similar  for  both  data  sets,  although  the  frequency  of  oscillations  is  greater  for  the  data  from  the  agar  target. 


102104  data.  TDR.  4096  ave. 


Figure  27:  TDR  measurements  with  and  without  an  agar  target  in  the  antenna.  The  blue  curve  in  the  lower  plot  is  the 
subtracted  signal  from  two  nonacoustic  TDR  data  sets,  taken  under  different  experimental  conditions.  The  green  curve  in  the 
lower  plot  is  the  subtracted  signal  from  acoustic  and  nonacoustic  data  sets  with  an  agar  target. 


We  have  previously  seen,  with  the  agar  target  experiments,  that  the  long  term  oscillations  in  the  subtracted 
data  remain  constant  whether  the  timing  circuit  is  externally  driving  the  signal  generator  or  the  signal 
generator  is  at  a  fixed  repetition  rate  (data  not  shown  in  this  report).  The  results  of  these  experiments  would 
seem  to  indicate  that  the  main  features  of  the  long-term  oscillations  in  the  subtracted  data  are  not  the  result 
of  acoustic  events,  but  an  artifact  of  the  transducer/timer  circuitry  interacting  with  the  electromagnetic 
signal  in  the  case  when  the  transducer  circuit  is  connected  to  the  timing  circuit. 
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7  COMPUTATIONAL  STUDIES 


As  the  above  experimental  results  indicate,  we  are  unable  to  detect  reflections  of  the  microwave  EM  pulse 
from  the  acoustic  wave  in  the  agar  material.  Previous  computational  work  illustrated  the  type  of  reflections 
that  would  be  observed  with  pressure  dependent  Debye  parameters  [1].  However,  the  parameter  values  used 
in  that  study  were  chosen  in  order  to  produce  a  measurable  effect  and  are  not  representative  of  a  particular 
material.  In  this  section  we  seek  to  obtain  an  estimate  for  the  pressure  dependence  of  the  dielectric  response 
of  agar  and  then  turn  again  to  computational  simulations,  this  time  to  attempt  to  explain  the  failure  to 
detect  EM  reflections  in  our  experimental  wave  guide. 


7.1  Dielectric  Properties  of  Agar 


Miura  et  al.  [13]  model  the  permittivity  of  agar  using  the  Havriliak-Negami  equation 

€sm  ^oo  ^sh  ^oo  /rj  -i  \ 

6  ~  600  =  (1  +  +  1  +  (iw^  ’  1  J 

where  u>  is  the  frequency  (Hz),  and  the  subscripts  m  and  h  refer  to  the  middle  and  high  frequency  regimes, 
respectively.  They  find  the  following  parameter  values  for  agar  (1%  by  weight)  esm  =  9.0,  Tm  =  1.82  x  10-8, 
am  =  0.94,  f3m  =  1.0,  esh  =  77.0,  Th  =  8.51  x  1CT12,  (3h  =  1.0,  and  =  5.0. 


Not  too  surprisingly,  the  permittivity  of  agar  closely  resembles,  in  the  high  frequency  limit,  that  of  water, 
which  has  a  single  molecular  relaxation  term  [8] 


£  ^OO 


^ s  ^oo 

1  +  (' iojT )  ’ 


(7.2) 


where  es  =  78.47  is  the  static  permittivity,  =  5.2  is  the  permittivity  in  the  high  frequency  limit,  and 
r  =  8.27  x  10  12  os  the  relaxation  time  for  molecular  orientation. 


7.2  Pressure  Dependence  of  the  Permittivity  of  Water 


Since  we  are  unable  to  find  data  in  the  literature  that  would  allow  us  to  estimate  the  pressure  dependence  of 
agar,  we  can  find  similar  values  for  water,  which  it  closely  resembles.  In  this  section  we  present  the  results 
of  a  literature  search  on  the  pressure  dependence  of  the  Debye  parameters  for  water. 

Following  the  work  of  Albanese  et  al.[  1]  we  first  approximate  the  pressure  dependence  of  the  Debye  parameters 
to  be  a  linear  perturbation  about  the  parameter  values  at  atmospheric  pressure 

t(p)  =  t0  +  tiTp 

es(p)  =  is,  o  +  KsP  (7.3) 

£oo  (p)  =  Coo.O  +  Koo p, 

where  p  is  the  pressure  and  kt,ks,  and  k oo  are  the  coefficients  of  pressure  for  r,  es,  and  e^,  respectively. 

We  were  unable  to  find  any  data  relevant  to  the  pressure  dependence  of  in  our  literature  search.  However, 
we  did  find  information  on  ks  and  kt,  which  we  present  below. 


7.3  Estimate  for  ks 


There  are  a  number  of  sources  in  the  literature  from  which  one  can  obtain  estimates  for  the  pressure 
dependence  of  es.  We  can  obtain  a  rough  estimate  for  (des/dP)\T  from  published  values  in  the  CRC 
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Handbook  ([8],  pg.  6-15).  Linearly  interpolating  from  data  at  0.1  and  1  MPa  and  295°K  and  300°K  ,  we  get 
room  temperature  T  =  298  °K  values  of  es  =  78.47  and  78.50  at  P  =  0.1  and  1.0  MPa,  respectively.  These 
es  values  provide  a  crude  estimate  of  ( des/dP)\T  =  ns  =  3.33  x  10-2  MPa-1  =  3.33  x  10-8  Pa-1  at  room 
temperature  and  atmospheric  pressure  ( P  =  0.101325  MPa). 

Fernandez  et  al.[ 9]  determined  that  es  =  78.41  and  (des/dP)\T  =  3.74  x  10-2  MPa-1  =  3.74  x  10-8  Pa-1 
at  T  =  298.144  °K  and  atmospheric  pressure. 

Floriano  and  Nascimento  [10]  extend  the  work  of  Bradley  and  Pitzer  [6]  to  more  accurately  fit  the  static 
dielectric  constant  es  to  published  data  [8]  at  low  temperatures  and  high  pressures  (0-2000  MPa).  Their 
model  for  the  static  dielectric  constant  at  room  temperature  (298°K  )  as  a  function  of  pressure  is  given  by 

es(P)  =  eR  +  a0ln  (  )  >  (7-4) 

\al  +  Pr  J 

where  Pr  is  a  reference  pressure,  in  their  case  10  MPa,  P  is  the  pressure  in  MPa,  and  eR  represents  the 
value  of  the  static  dielectric  constant  at  the  reference  pressure.  At  T  =  298°K  ,  ao  =  14.1113,  ai  =  341.5902 
MPa,  and  eR  =  78.85. 


Since  we  are  interested  in  a  linearized  model  (7.3),  we  can  linearize  (7.4)  about  Pq  =  0.1  MPa  by  writing 
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-  a  ,  P  (7-10) 
ai  +  Jo 

1  atmosphere  =  0.101  MPa,  gives  es  =  78.44  and  ns  =  4.13  x  10-2  MPa-1  =  4.13  x  10-8 


7.4  Estimate  for  kt 

Using  data,  Shcherbakov  [16]  estimates  the  pressure  dependence  of  r  to  be 

t(P)  =  T0  -  k*P*  (s), 


where  P*  is  the  pressure  in  MPa, 

*  f  3-48  x  lO-12??*  ^ 

^  =  l - T - ) 

and  ij*  is  the  viscosity  in  centipoise. 

The  above  equations  can  be  rewritten  as 


=  1.05  x  10“14 


(s/MPa), 


where  P  is  the  pressure  in  Pa  and 


t(P)  =  tq  —  ktP  (s), 


kt  =  1.05  x  10”20  (s/Pa). 


(711) 


(7.12) 


(7.13) 

(7.14) 
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7.5  Summary  of  Literature  Search 


In  summary,  we  have  obtained  estimates  for  the  pressure  dependence  of  the  Debye  parameters  es  and  r  for 
water  from  data  available  in  the  literature.  No  data  was  found  for  the  pressure  dependence  of  6^.  Based  upon 
the  measured  frequency  dependence  of  the  permittivity  of  agar  [13]  we  expect  that  the  pressure-dependence 
of  the  Debye  parameters  for  water  would  be  a  good  first  estimate  for  the  pressure-dependent  behavior  of 
the  permittivity  of  agar.  Values  for  the  pressure  coefficients  (kt  =  —1.05  x  10-2°  and  ks  =  3.74  x  10-8) 
that  were  obtained  from  the  literature  are  much  smaller  than  the  values  used  in  previous  simulations  [1]. 
Therefore,  we  suggest  that  our  experimental  conditions  are  probably  not  sufficient  to  see  reflections  of  the 
EM  pulse  from  the  pressure  wave.  (We  estimate  the  magnitude  of  the  pressure  waves  from  the  transducer 
to  be  on  the  order  of  0.1  Pa.) 

The  table  below  summarizes  the  Debye  parameter  values  for  water  at  atmospheric  pressure  and  estimates 
from  the  literature  for  the  (linear)  coefficients  of  pressure. 


Table  1:  Parameter  values  for  water  at  room  temperature  and  atmospheric  pressure 


to  [s] 

e«,o 

^oo,0 

kt  [s-Pa  x] 

ks  [Pa-1] 

CRC  Handbook 

8.27x  10-12 

78.47 

5.2 

3.3xl0-8 

Fernandez  et  al. 

78.41 

3.74xl0-8 

Floriano  and  Nascimento 

78.44 

4.13xl0-8 

Shcherbakov 

8.2xl0-12 

-1.05xl0-20 

7.6  Simulations 


We  present  numerical  simulations  for  a  Debye  medium  that  is  similar  to  water.  For  signal  bandwidths  in 
the  microwave  regime  the  dispersive  properties  of  water  are  usually  modeled  by  a  Debye  equation  having  a 
single  molecular  relaxation  term.  The  mean  values  eSio>  eoo,o>  and  r0  for  this  Debye  model  are  given  by 

eSjo  =  78.44  (relative  static  permittivity), 

Coo,0  =  5.2  (relative  high  frequency  permittivity), 

t0  =  8.2  x  10-12  seconds,  ’  (7-15) 

a  =  1  x  10-5  mhos/meter, 

where  a  is  the  conductivity  of  water. 

Following  the  notation  in  [3]  we  represent  the  pressure  dependent  parameters  es,  Cqo,  and  t  as  a  mean  value 
plus  a  perturbation  that  is  proportional  to  the  pressure,  with  the  coefficients  of  pressure  designated  as  ks, 
Kqo  and  kt,  respectively.  Based  upon  the  above  literature  search  and  analysis,  we  estimate  the  following 
values  for  the  coefficients  of  pressure  dependence:  ks  =  3.74  x  10  8Pa  1  and  kt  =  —1.05  x  10  20  s/Pa,  with 
Koo  undetermined.  As  demonstrated  in  [3],  the  acoustic  reflection  is  not  sensitive  to  changes  in  either  Koo 
or  kt,  but  is  sensitive  to  changes  in  ks.  Based  on  these  observations  and  the  low  value  of  kt,  we  chose  the 
following  values  for  our  simulations 

ks  =  3.74  x  10-8  Pa-1, 

Koo  =0.0,  ,  (7.16) 

kt  =0  s/Pa. 
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PEC  wall 


The  computational  domain  is  similar  to  that  of  [3]  and  is  defined  as  follows.  We  take  Xq  =  (0,0.1), 
Za  =  (0,0.15)  and  Zn  =  (0.15,0.2).  The  number  of  nodes  along  the  2-axis  is  taken  to  be  320  and  the 
number  of  nodes  along  the  x-axis  is  taken  to  be  160.  The  spatial  step  size  in  both  the  x  and  2  directions  is 
Ax  =  A 2  =  h  =  0.1/160.  From  the  CFL  condition  with  the  Courant  number  rj cn  =  1/2  we  obtain  the  time 
increment  to  be  At  «  1.0417  pico  seconds.  The  central  frequency  of  the  input  source  as  described  in  (7.19) 
is  3.0  GHz  and  based  on  the  speed  of  light  in  air,  cq  =  3  x  108  m/s,  we  calculate  the  corresponding  central 
wavelength  to  be  Ac  =  (27tco)/w  =  0.1  meters.  The  source/observation  antenna  is  half  a  wavelength  long 
and  is  placed  at  (. xi,X2 )  x  zc,  with  2C  =  0,  xi  =  0.025  and  X2  =  0.075.  We  use  PML  layers  that  are  half  a 
wavelength  thick  on  all  fours  sides  of  the  computational  domain  as  shown  in  Figure  28.  The  reflections  of 
the  electromagnetic  pulse  at  the  air-Debye  interface  and  from  the  acoustic  pressure  wave  are  recorded  at  the 
center  of  the  antenna  (xc,zc),  with  xc  =  0.05,  at  every  time  step.  The  component  of  the  electric  field  that 
is  of  interest  here  is  the  Ex  component  of  E.  Thus  our  “observations”  for  the  reflected  signals  are 
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xc ,  zc ,  q  ) } 

(5k  5k  sk  sk  5k 

es,0>eoo,0>r0^ 
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n— 1 

*£o> 
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O 


(7.17) 


where  E  is  the  computed  solution  to  Maxwell’s  equation  as  given  in  [3].  In  particular,  we  use  modified 
Maxwell’s  equations  which  govern  the  electric  field  E  and  the  magnetic  field  H  in  a  domain  O  with  charge 
density  p.  Thus  we  first  consider  the  system 


(i)  ^  +  J-VxH  =  0,  in  (0,  T )  x  fi, 

(ii)  —  +  V  x  E  =  0,  in  (0,  T)  x  Q, 
at 

(iii)  V  ■  D  =  p,  in  (0,  T)  x  f l, 

(iv)  V  •  B  =  0,  in  (0,  T)  x 

(v)  E  x  n  =  0,  on  (0,  T)  x  dfl, 

,  (vi)  E(0,x)  =  0,  H(0,x)  =  0,  in  O. 


We  have  J  =  Jc  +  Js,  where  Jc  is  a  conduction  current  density  and  Js  is  the  source  current  density.  We 
assume  only  free  space  (actually  the  antenna  in  our  example)  can  have  a  source  current,  and  Jc  is  only  found 
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in  the  dielectric  material.  The  electromagnetic  input  source  Js  has  the  form 


J s{t,  x,  z)  =  I(Xl>X2)S(z)  sin (vc(t  -  3 10))  exp 

1 


t  —  3to 


to 


to  — 


sec,  tuc  =  6ir  x  109  rad/sec,  fc  =  3.0  x  109  Hz. 


2tt  x  109 

The  Fourier  spectrum  of  this  pulse  has  even  symmetry  about  3.0  GHz. 


(7.19) 


To  the  system  (7.18),  we  add  variables  to  account  for  the  PML  layers  surrounding  the  computational  domain- 
see  [3]  for  complete  details.  Within  the  dielectric  medium  we  have  constitutive  relations  that  relate  the  flux 
densities  D,  B  to  the  electric  and  magnetic  fields,  respectively.  We  have 

(i)  D  =£0E  +  P/d, 

(ii)  B  =  ^0H,  (7.20) 

,  (in)  J c  =  ctE/d, 

In  (7.20),  Id  denotes  the  indicator  function  on  the  Debye  medium.  Thus,  Jc  =  0  in  air.  The  electric 
polarization  P  is  given  by 

P  =  Pi  +  PR  =  e0\-E  +  PR,  (7.21) 

and  hence  the  constitutive  law  (7.20,  i)  becomes 

D  =  e0erE  +  PR,  (7.22) 

where  er  =  (1  +  x)  is  the  relative  permittivity  of  the  dielectric  medium.  The  Debye  model  is  represented  in 
differential  form  as 


tPr  +  PR  =  £o(es  -  e00)E, 


D  —  eoCooE  +  PR, 


(7.23) 


inside  the  dielectric,  whereas  PR  =  0,  =  1  in  air,  and  =  er  inside  the  dielectric.  We  will  henceforth 

denote  PR  by  P.  In  equation  (7.23),  the  parameters  es  and  denote  the  static  relative  permittivity,  and 
the  value  of  permittivity  for  an  extremely  high(ss  oo)  frequency  field,  respectively. 

The  pressure  p  will  be  assumed  to  have  the  form 

P{t,z)  =  I(Z2,Z2+xP)\p\sm  (u)p[t+  -]^,  (7.24) 

where  z^  £  Z R  with  z\  <  z^-  The  terms  wp,  Xp  \p\  and  cp  denote  the  acoustic  frequency,  wavelength, 
pressure  magnitude  and  speed  respectively.  The  windowed  acoustic  pressure  wave  as  defined  in  (7.24)  has 
the  parameter  values  top  =  6.07T  x  105  rad/sec,  cp  =  1500  m/s,  and  thus,  Xp  =  0.005.  Also  \p\  =  1  Pa. 
The  location  of  the  pressure  region  is  in  the  interval  (22,^2  +  Xp)  =  (0.175,0.18).  Further  details  on  these 
simulations  can  be  found  in  [3]. 

Since  only  na\p\  enters  into  these  calculations,  it  is  only  the  product  that  matters.  In  Figure  29  we  plot  the 
ratio  of  the  maximum  amplitude  of  the  acoustic  reflection  to  the  maximum  amplitude  of  the  input  source 
for  different  values  of  the  product  /cs|p|.  The  other  pressure  coefficients,  Koo  and  nT  are  set  to  be  zero.  From 
our  simulations  we  see  that  this  ratio  |  Ea  \  / 1  Einput  is  independent  of  the  maximum  amplitude  of  the  source 
term,  i.e.,  \Ea\  scales  according  to  the  value  of  |Einput|.  We  can  see  in  Figure  29  that  the  ratio  | Ea \ / 1 Einput 
is  linearly  dependent  on  the  ratio  ns\p\  for  values  larger  than  10.  We  note  that  E  =  E{ns\p\)  is  not  linear  in 
Ks\p\  even  though  we  make  the  linear  approximations  of  (7.3)  for  r,  es  and  as  a  function  of  p. 
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Figure  29:  Plot  of  the  ratio  l-EUl/l-Einputl  versus  the  product  Ks|p|. 


We  plot  the  ratio  of  the  maximum  amplitude  of  the  reflected  wave  to  the  maximum  amplitude  of  the  input 
wave,  since  the  amplitude  of  the  reflected  wave  scales  with  the  amplitude  of  the  input  wave.  We  plot  these 
values  as  a  function  of  the  product  /cs|p|,  since  the  strength  of  the  reflection  scales  with  this  product.  In 
this  way,  the  curve  can  be  used  in  a  universal  manner.  Knowing  the  pressure  dependence  ks  of  a  material, 
one  can  determine  the  magnitude  of  the  pressure  wave  necessary  to  achieve  a  desired  signal-to-input  ratio. 
Similarly,  if  one  knows  the  inherent  noise  in  making  observations  of  the  reflected  waves,  one  can  determine 
the  size  of  the  input  source  needed  to  achieve  a  satisfactory  signal-to-noise  ratio,  i.e.,  to  produce  a  detectable 
effect. 

Figure  30  depicts  the  same  data  as  Fig.  29,  except  on  a  log  x  scale.  With  this  view,  it  can  be  seen  that  the 
values  for  the  ratio  |.EU|/|-E'input  plateau  for  small  values  of  Ks\p\.  Simulations  with  different  grid  refinements 
show  that  this  plateau  tends  toward  zero  as  the  mesh  size  decreases. 

We  can  use  these  curves  to  estimate  the  amplitude  of  the  pressure  wave  that  we  would  need  in  order  to 
observe  a  reflection  of  the  EM  pulse  in  our  experiments.  The  magnitude  of  the  incident  electric  field  in  the 
experiment  is  roughly  15  V /m.  If  we  look  for  the  amplitude  of  the  reflected  pulse  to  be  roughly  15  mV /m 
(in  order  to  be  above  the  noise),  then  we  would  look  for  a  ratio  |Fbi|/|.EinpUt|  =  1  x  10-3.  This  would 
require  that  ks\p\  ~  30.  Since  ns  =  3.74  x  10-8  for  water  (agar),  we  estimate  that  we  would  need  a  pressure 
amplitude  \p\  «  8  x  108  Pa,  or  8000  atmospheres. 
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Figure  30:  Plot  of  the  ratio  |i?A|/|F^input|  versus  the  product  log(/cs|p|). 


8  CONCLUSION 


In  summary,  we  have  constructed  a  workable  experimental  antenna  where  we  can  launch  EM  waves  at  an 
agar  target  and  make  TDR  measurements.  The  circuitry  and  timing  necessary  to  initiate  an  acoustic  wave 
and  synchronize  it  with  the  EM  pulse  have  been  built.  We  have  carried  out  TDR  experiments  with  different 
delay  times,  pulse  widths,  and  target  thicknesses. 

To  date,  there  is  no  evidence  of  EM  reflection  from  the  acoustic  wave  front  for  the  range  of  inputs  in  which 
our  experimental  antenna  can  be  operated.  Experimental  results  are  complicated  due  to  jitter  in  the  time 
base  which  results  in  slight  offsets  and  spikes  in  the  subtracted  signals.  These  spikes  in  the  subtracted  signal 
are  also  observed  when  the  input  pulse  reflects  from  the  feed  of  the  antenna,  further  obscuring  results  in  the 
vicinity  of  interest.  In  addition,  the  long-term  oscillations  that  have  been  observed  in  the  subtracted  data 
seem  not  to  be  the  result  of  acoustic  events,  but  an  artifact  of  the  transducer/timer  circuitry  interacting 
with  the  electromagnetic  signal  in  the  case  when  the  transducer  circuit  is  connected  to  the  timing  circuit. 

Evaluation  of  the  expected  pressure  dependence  of  the  agar  indicates  that  it  would  be  difficult  to  observe  EM 
reflection  from  the  acoustic  wave  front  with  our  current  transducer,  where  the  magnitude  of  the  pressure  wave 
is  estimated  to  be  on  the  order  of  0.1  Pa.  range.  A  suitable  source  for  acoustic  waves  must  satisfy  a  difficult 
set  of  constraints.  It  must  be  much  more  powerful  than  our  current  transducer,  must  be  remotely  controlled 
for  synchronization  of  the  EM  and  acoustic  signals,  must  produce  a  high  frequency  pulse  that  produces  a 
sharp  wave  front  in  the  material,  and,  ideally,  would  not  interact  with  the  EM  field  in  the  antenna  (no  wires 
or  metal  parts). 
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